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Abstract 

We show how to express the information contained in a Quantum Bayesian (QB) net 
as a product of unitary matrices. If each of these unitary matrices is expressed as 
a sequence of elementary operations (operations such as controUed-nots and qubit 
rotations), then the result is a sequence of operations that can be used to run a 
quantum computer. QB nets have been run entirely on a classical computer, but one 
expects them to run faster on a quantum computer. 
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Introduction 



Quantum Bayesian (QB) Nets|I[]-0 are a method of modelling quantum sys- 
tems graphically in terms of network diagrams. In this paper, we show how to express 
the information contained in a QB net as a product of unitary operators. In another 
paper 1^, we have presented a method for reducing any unitary operator to a sequence 
of elementary operations (SEO). Such SEOs can be used to manipulate an array of 
quantum bits, a quantum computer. [^]-|]^ Thus, combining the results of this pa- 
per and those of Ref. IQ, one can reduce a QB net into a SEO that can be used to 
run a quantum computer. Of course, QB nets can and have been run entirely on a 
classical computer 0. However, because of the higher speeds promised by quantum 
parallelism, one expects QB nets to run much faster on a quantum computer. 

The process of reducing a QB net to a SEO is analogous to the process of "com- 
piling source code" for classical computers. With classical computers, one writes a 
computer program in a high level language (like Fortran, C or C-I-+). A compiler 
then expresses this as a SEO for manipulating bits. In the case of quantum com- 
puters, a QB net may be thought of as a program in a high level language. In the 
future, programs called "quantum compilers" will be widely available. These will run 
(initially) on classical computers. They will take whatever QB net we give it, and 
re-express it as a SEO that will then be used to control an array of quantum bits. 
These quantum compilers will take the QB net entered by the user and automatically 
add to it quantum error correction code[0 and various speed optimizations. 

QB nets are to quantum physics what Classical Bayesian (CB) nets]^ are to 
classical physics. CB nets have been used very successfully in the field of artificial 
intelligence (AI). In fact, even Microsoft's Office Suite contains CB net code[0. Thus, 
we hope and expect that some day QB nets, running on quantum computers, will 
be used for AI applications. In fact, we believe that quantum computers are ideally 
suited for such applications. First, because AI tasks often require tremendous power, 
and quantum computers seem to promise this. Second, because quantum computers 
are plagued by quantum noise, which makes their coherence times short. There 
are palliatives to this, such as quantum error correction[0. But such palliatives 
come at a price: a large increase in the number of steps. The current literature 
often mentions factoring a large number into primes as a future use of quantum 
computers. However, due to noise, quantum computers may ultimately prove to be 
impractical for doing long precise calculations such as this. On the other hand, short 
coherence times appear to be a less serious problem for the types of calculations 
involved in AI. The human brain has coherence times too short to factor a 100 digit 
number into primes, and yet long enough to conceive the frescoes in the Sistine 
Chapel. We do not mean to imply that the human brain is a quantum computer. 
An airplane is not a bird, but it makes a good flyer. Perhaps a quantum computer, 
although not a human brain, can make a good thinker. 
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Review 



We begin by presenting a brief review of QB nets. For more information, see 
Ref.i-i. 

In what follows, we use the following notation. We define Z^f, = {a, a + 
1, . . . ,b} for any integers a and b. S{x, y) equals one if x = y and zero otherwise. For 
any finite set S, \S\ denotes the number of elements in S. 

We call a graph (or a diagram ) a collection of nodes with arrows connecting 
some pairs of these nodes. The arrows of the graph must satisfy certain constraints. 
We call a labelled graph a graph whose nodes are labelled. A QB net consists of 
two parts: a labelled graph with each node labelled by a random variable, and a 
collection of node matrices, one matrix for each node. These two parts must satisfy 
certain constraints. 

An internal arrow is an arrow that has a starting (source) node and a different 
ending (destination) one. We will use only internal arrows. We define two types of 
nodes: an internal node is a node that has one or more internal arrows leaving it, and 
an external node is a node that has no internal arrows leaving it. It is also common 
to use the terms root node or prior probability node for a node which has no incoming 
arrows (if any arrows touch it, they are outgoing ones). 

We restrict our attention to acyclic graphs; that is, graphs that do not contain 
cycles. (A cycle is a closed path of arrows with the arrows all pointing in the same 
sense.) 

We assign a random variable to each node of the QB net. (Henceforth, we 
will underline random variables. For example, we might write P{x = x) for the 
probability that the random variable x assumes the particular value x.) Suppose the 
random variables assigned to the nodes For each j G Zi^^, the 

random variable Xj will be assumed to take on values within a finite set Ej called the 
set of possible states ofxj. 

If S* = {ki, k2, ■ ■ ■ , k\s\} C Zi^N, and ki < k2 < ■ ■ ■ < k\s\, define {x-)s = 
(xfc^, Xfcj, ■ ■ ■ , Xfc|g|) and {x.)s = {xj^j^jX/.^, - ' ' j^k^s\)- Sometimes, we also abbreviate 
{x-)ziN the vector that includes all the possible Xj components) by just x., and 
(^•)^i,iv by just X. . 

Let Zext be the set of all j G Zi^n such that Xj is an external node, and let 
Zint be the set of all j G tv such that Xj is an internal node. Clearly, Z^xt and Zint 
are disjoint and their union is Zi^j^. 

Each possible value x. of x. defines a different net story. For any net story x-, 
we call {x-)zi„t the internal state of the story and {x-)z^^t its external state. 

For each net story, we may assign an amplitude to each node. Define Sj to be 
the set of all k such that an arrow labelled Xk (i.e., an arrow whose source node is 
Xj^) enters node Xj. We assign a complex number v4j[2:j|(z.)sj to node Xj. We call 
ylj[a;j|(x.)5j the amplitude of node Xj within net story x-. 



3 



The amplitude of net story x., call it A{x.), is defined to be the product of all 
the node amplitudes 74j[xj|(x.)5j for j G Zi j^. Thus, 

A{x.)= n M^j\i^-)s,]- (1) 

The function Aj with values (a;.)^^] determines a matrix that we call the 

node matrix of node Xj. xj is the matrix's row index and {x-)sj is its column index. 

Method 



2 qbits 



1 qbits 



2 qbits 




1 qbits 



1 qbits 



Fig.l QB net for Teleportation. This figure also shows the number of quantum or 

classical bits carried by each arrow. 

One can translate a QB net into a SEO by performing the following 3 steps: 
(1) Find eras, (2) Insert delta functions, (3) Find unitary extensions of era matrices. 
Next we will discuss these steps in detail. We will illustrate our discussion by using 
Teleportation as an example. Figure 1 shows a QB net for Teleporation. This 
net is discussed in Ref. 0. Reference |TT| gives a SEO, expressed graphically as a 
qubit circuit, for Teleportation. It appears that the author of Ref.|[ll|] obtained his 
circuit mostly by hand, based on information very similar to that contained in a QB 
net. This paper gives a general method whereby such circuits can be obtained from 
a QB net in a completely mechanical way by means of a classical computer. 
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Step 1: Find eras 




Fig. 2 Root node eras for Teleportation net. 



The root node eras of a graph are defined as follows. Call the original graph 
Graph(l). The first era Ti is defined as the set of all root nodes of Graph(l). Call 
Graph(2) the graph obtained by erasing from Graph(l) all the Ti nodes and any 
arrows connected to these nodes. Then T2 is defined as the set of all root nodes 
of Graph(2). One can continue this process until one defines an era T]-r| such that 
Graph(|T| + 1) is empty. (One can show that if Graph(l) is acyclic, then one always 
arrives at a Graph(|T| + l) that is empty.) For example, Fig. 2 shows the root node eras 
for the Teleportation net Fig.l. Let T represent the set of eras: T = {Ti, T2, ■ ■ ■ , 
Note that Ta C Zi^n for all a G Zi^\q-\ and the union of all Ta equals Zi^n. In 
mathematical parlance, the collection of eras is a partition of Zi^n- 
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Fig. 3 External node eras for Teleportation net. 



Rather than defining eras by (1) removing successive layers of root nodes, one 
can also define them by (2) removing successive layers of external nodes. We call this 
second type of era, the external node eras of the graph. For example. Fig. 3 shows the 
external node eras of the Teleportation net Fig.l. 

This process whereby one classifies the nodes of an acyclic graph into eras is a 
well know technique referred to as a chronological or topological sort in the computer 
literature [p^] . 

Henceforth, for the sake of definiteness, we will speak only of root node eras. 
The case of external node eras can be treated similarly. 

Suppose that a G Zi^\t\. The arrows exiting the a'th era are labelled by {x.)t^. 
Those entering it are labelled by {x.)y^, where Fq is defined by F^ = [jjeTa '^j- Note 
that the a'th era node is only entered by arrows from nodes that belong to previous 
(not subsequent) eras so F^ C Ta-i U . . . U T2 U Ti. The amplitude Ba of the a'th era 
is defined as 



i?a[(x.)Tj(a;.)rJ = U A,[x,\{x.)s,] . 



(2) 



jeTa 

The amplitude A{x.) of story x- is given by 



A{x.) = l[B, 



(3) 



a=l 
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For example, for Teleportation we get from Fig. 2 

Bi{xi,X4) = Ai{xi)A4{x4) , 



(4a) 



B2{x2,X3\xi) = A2{x2\xi)A3{x3\xi) , (4b) 

B3{X5\X2,X4) = A5{X5\X2,X4) , (4c) 
Bi{xG\x3,X5) = Ae{xG\x3,X5) , (4d) 

and 

A{x.) = B4B3B2B, . (5) 
Step 2: Insert delta functions 

The Feynman Integral FI for a QB net is defined by 

FI[{x-)z.J^ E M^-)- (6) 

Note that we are summing over all stories x. that have {x.)z^^t their external state. 
We want to express the right side of Eq. (6) as a product of matrices. 

Consider how to do this for Teleportation. In that case one has 

F/(x6)= E B4BSB2B,, (7) 

where the Ba are given by Eqs(4). The right side of Eq.(7) is not ready to be expressed 
as a product of matrices because the column indices of Ba+i and the row indices of 
Ba are not the same for all a e Zi |r|-i. Furthermore, the variable X3 occurs in B4 
and B2 but not in B3. Likewise, the variable X4 occurs in B3 and Bi but not in B2. 
Suppose we define B^ for a e Zi^\t\ by 

Bi{xl,xl) = Bi{xl,xl) , (8a) 

B2{x2, x1, x1\x\, x\) = B2{x2, x1\x\)5{x1, x\) , (8b) 

B3{xl,xl\xl,xl,x1) = B3{xl\xl,x1)S{xl,xl) , (8c) 

B4{xq\xI xl) = B4{xq\xI xl) . (8d) 

Then 
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1 ) 



(9) 



interm 



where we sum over all intermediate indices; i.e., all except x^. Contrary to Eq.(7), 
the right side of Eq.(9) can be expressed immediately as a product of matrices since 
now Ba+i column indices and Ba row indices are the same. The purpose of inserting 
a delta function of into B^ is to allow the system to "remember" the value of x^ 
between non-consecutive eras T4 and T2. Inserting a delta function of X4 into B2 
serves a similar purpose. 



1 ^ 









X 



Fig.4 Example of a QB net in which an external node is not in the final era. 

In the Teleporation net of Fig.l, the last era contains all the external nodes. 
However, for some QB nets like the one in Fig.4, this is not the case. For the net of 
Fig.4, 



Bi{xi) = Ai{xi) , 

B2{x2,X3\xi) = As{x3\xi)A2{x2\Xi] 



'10a) 



(10b) 
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(10c) 



54(^5 1 0:4) = A (3:5 1 2^4) • (lOd) 

Even though node X2 is external, the variable X2 does not appear as a row index in 
S4. Suppose we set 





Blixl) 


= B^{x\) , 


(11a) 




B2{xl,xl\x[) 


= B2{xl,xl\x\) , 


(lib) 




U/2) •''4|'t'2) ■''3) / 


= B3{xl\xl)6{xl, xl) , 


(11c) 


S4 


(X2, ^5 |X2, 3^4) = 


= 54(2:5 a;4)(5(x2, 0:2) . 


(lid) 



Then 

FI{X2,X5)^ J2 BAB3B2B1, (12) 

interm 

where we sum over all intermediate indices; i.e., all x| except X2 and X5. Contrary to 
i?4, the rows of i?4 are labelled by the indices of both external nodes X2 ^-^d x^. 

This technique of inserting delta functions can be generalized as follows to 
deal with arbitrary QB nets. For j e ^i,Ar, let aminU) be the smallest a e -Z^ijri such 
that Xj appears in Hence, a^ini,]) is the first era in which Xj appears. If is 
an internal node, let amax{j) be the largest a such that Xj appears in Ba (i.e., the 
last era in which Xj appears). If Xj is an external node, let amaxU) = |^| + 1- ^oi 
a e Zi^\r\, let 

A„ = {j e ^i,jv|amm(j) <a< amaxij)} , (13) 

Ba = 5j(x?)Tj(x?-^)rJ n ■ (14) 

In Eq.(14), should be identified with Xj and x° with no variable at all. Equation(6) 
for FI can be written in terms of the Ba functions: 

FI[{x.)z^J= E B^r\...B2B,, (15) 

interm 

where the sum is over all intermediate indices (i.e., all Xj for which a ^ For all 

a, define matrix Ma so that the x, y entry of Ma is Ba{x\y). Define M to be a column 
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vector whose components are the values of FI for each external state. Then Eq.(15) 
can be expressed as: 

M = M\r\ . . . M2M1 . (16) 

The rows of the column vector M are labelled by the possible values of {x.)z^^f The 
rows of the column vector Mi are labelled by the possible values of where Ti 

is the set of root nodes. 

Step 3: Find unitary extensions of era matrices 

So far, we have succeeded in expressing FI as a product of matrices M^, 
but these matrices are not necessarily unitary. In this step, we will show how to 
extend each Ma matrix (by adding rows and columns) into a unitary matrix Ua- The 
techniques of Ref.[^ will then be applicable to each matrix Ua- 

By combining adjacent M^'s, one can produce a new, smaller set of matrices 
Ma- Suppose the union of two consecutive eras is also defined to be an era. Then 
combining adjacent M^'s is equivalent to combining two consecutive eras to produce a 
new, smaller set of eras. We define a breakpoint as any position a G between 
two adjacent matrices M^+i and Ma- Combining two adjacent M^'s eliminates a 
breakpoint. Breakpoints are only necessary at positions where internal measurements 
are made. For example, in Teleportation experiments, one measures node X3, which 
is in era T3. Hence, a breakpoint between M4 and M3 is necessary. If that is the 
only internal measurement to be made, all other breakpoints can be dispensed with. 
Then we will have M = M^M[ where M^ = M4, M[ = M3M2M1. If no internal 
measurements are made, then we can combine all matrices Ma into a single one, and 
eliminate all breakpoints. 

We will henceforth assume that for all a G Zi^\q-\, the columns of Ma are 
orthonormal. If for some Oq ^ Miq does not satify this condition, it may be 

possible to "repair" Mag so that it does. First: If a row (3 of Ma^^i is zero, then 
eliminate the column f3 of Mq^, and the row (3 of Ma^-i- Next: If a row (3 of the 
column vector Ma^-i - - - M2M1 is zero, then flag the column /3 of Ma^- The flagged 
columns of M^y can be changed without affecting the value of M. If the non-fiagged 
columns of Ma^ are orthonormal, and the number of columns in Mag does not exceed 
the number of rows, then the Gram Schmidt method, to be discussed later, can be 
used to replace the flagged columns by new columns such that all the columns of the 
new matrix Mag are orthonormal. If it is not possible to repair Mag in any of the 
above ways (or in some other way that might become clear once we program this), 
one can always remove the breakpoint between Ma^+i and Mag- 
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V = U3 

Fig.5 Dimensions of matrices Ma and of their unitary extensions 

Consider Fig.5. We will call da the number of rows of matrix Ma and da-i its 
number of columns. We define D and Ns by 

L> = max{4|l < a < |T|} (17) 

Ns = min{2^|i G Zi,oo, D < 2'} . (18) 

Let da = Ns — da for all a. For each a 7^ 1, we define [/„ to be the matrix that one 
obtains by extending Ma as follows. We append an da x da^i block of zeros beneath 
Ma and an A^^ x da-i block of gray entries to the right of Ma. By gray entries we 
mean entries whose value is yet to be specified. When a = 1, Mi be can extended in 
two ways. One can append a column of di zeros beneath it and call the resulting A^^ 
dimensional column vector vi. Alterntatively, one can append a column of di zeros 
beneath Mi and an Ns x (A^^ — 1) block of gray entries to the right of Mi, and call 
the resulting A^^ x A^^ matrix Ui. In this second case, one must also insert ei to 
the right of Ui. By ci we mean the Ns dimensional column vector whose first entry 
equals one and all others equal zero. Which extension of Mi is used, whether the 
one that requires ei or the one that doesn't, should be left as a choice of the user. 
Henceforth, for the sake of definiteness, we will will assume that the user has chosen 
the extension without the Ci. The other case can be treated similarly. Equation(16) 
then becomes 

v = Uir\... U3U2V1 , (19) 
where v is just the column vector M with d\r\ zeros attached to the end. Note that 



Ma... M2M1 





(20) 
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for all a G -Z'l |7-|, where the zero indicates a column of da zeros. 

To determine suitable values for the gray entries of the Ua matrices, one can 
use the Gram-Schmidt (G.S.) method [0. This method takes as input an ordered 
set 5* = {vi,V2, ■ ■ ■ ,vn) of vectors, not necessarily independent ones. It yields as 
output another ordered set of vectors S' = {ui,U2, . . . ,un), such that S' spans the 
same vector space as S. Some vectors in S' may be zero. Those vectors of 5" which 
aren't zero will be orthonormal. For r G 2'i tv, if the first r vectors of S are already 
orthonormal, then the first r vectors of S' will be the same as the first r vectors of 
5". Let Cj for j G Zi^jy^ be the j'th standard unit vector (i.e., the vector whose j'th 
entry is one and all other entries are zero). For each a G Zi^\q-\, to determine the gray 
entries of Ua one can use the G.S. method on the set S consisting of the non-gray 
columns of Ua together with the vectors ei,e2, ■ ■ - Cns- 
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FIGURE CAPTIONS: 

FiG.l QB net for Teleportation. This figure also shows the number of quantum or 
classical bits carried by each arrow. 

Fig. 2 Root node eras for Teleportation net. 

Fig. 3 External node eras for Teleportation net. 

Fig. 4 Example of a QB net in which an external node is not in the final era. 
Fig. 5 Dimensions of matrices and of their unitary extensions Ua- 



14 



